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Abstract 

We investigate the kinematics of deformations in two and three dimensional media by explicitly 
solving (analytically) the evolution equations (Raychaudhuri equations) for the expansion, shear 
and rotation associated with the deformations. The analytical solutions allow us to study the 
dependence of the kinematical quantities on initial conditions. In particular, we are able to identify 
regions of the space of initial conditions that lead to a singularity in finite time. Some generic 
features of the deformations are also discussed in detail. We conclude by indicating the feasibility 
and utility of a similar exercise for fluid and geodesic flows in flat and curved spacetimes. 
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I. INTRODUCTION 



Deformations (elastic or otherwise) under applied external stress is a well-studied subject 
It is known that the kinematics of such deformations can be understood by analysing 
the behaviour of the kinematical quantities, namely the isotropic expansion, the shear and 
the rotation or vorticity. However, as far as we know, most of these investigations do not 
explicitly address the question of evolution of the expansion, rotation or shear in time by 
solving the associated initial value problem for the system of differential equations involving 
these variables. We believe that such a study could shed light on the role of initial conditions 
(more precisely, the initial values of the expansion, rotation and shear at some time) in the 
kinematics of deformations. Additionally, many subtleties (integrals of motion, critical val- 
ues, appearance of singularities etc.) in the nature of the evolution of these quantities may 
emerge if one is able to obtain exact solutions or even scan the full domain of numerical so- 
lutions of the relevant coupled, nonlinear system of evolution equations involving expansion, 
shear and rotation. 

The equations we will discuss in this article are, essentially, the Raychaudhuri equations 
so prominent and well-known in the study of spacetime singularities in gravitation and 
cosmology In the context of gravitation, the central result that emerges from the 

Raychaudhuri equations, is the notion of geodesic focusing. The analysis which is based on 
an inequality (in the first order approach) or the theory of second order ordinary differential 
equations (in the second order approach) proves the seemingly simple statement: if gravity 
is attractive then geodesies must focus to a singularity in the congruence (which, may or may 
not be a genuine spacetime singularity). In the broader context of Riemannian geometry, 
focusing in a geodesic flow arises provided one obeys the appropriate convergence condition. 

In this article, we, however, will restrict ourselves to a background flat space (not space- 
time) and first deal with two spatial dimensions, and then with three dimensions. The 
family of worldlines (the so-called geodesic congruence) will be replaced, in our case (i.e. for 
deformable media), by the deformation vectors, with the time derivative of the deformation 
vector replacing the tangent vector to the geodesic flow. The geodesic equation used to gen- 
erate the geodesic flow, which is normally used as an input in deriving the Raychaudhuri 
equations (one can also write equations for non-geodesic flows too, but we do not consider 
them here) will now be replaced by a second order (in time) equation for the deformation 
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vector. In other words, we shall exactly follow the scheme of derivation of the Raychaud- 
huri equations in the context of General Relativity, keeping in mind the abovementioned 
parallels and restrictions. We mention here that by virtue of being a geometric statement, 
the Raychaudhuri equations are applicable to a myriad of diverse situations ranging from 
deformations in media, fluid flow, geodesic congruences etc. Our aim, as mentioned before, 
is to analyse, in all possible detail, the kinematics of deformations of a medium without and 
with stiffness and viscosity. In particular, one of the features we investigate is: how does 
the appearance of a singularity depend on the initial conditions on the expansion, shear or 
rotation. The paper is broadly divided into two parts. In Section II we talk about two di- 
mensional media, while in Section III we focus attention on three dimensional media. In each 
part, we write down the relevant equations involving the expansion, rotation and shear and 
find analytic solutions for these variables. Subsequently, we try to understand the solutions 
obtained by looking for generic features and their dependence on initial conditions. 

II. DEFORMATIONS IN TWO DIMENSIONAL MEDIA 
A. Expansion, rotation and shear 

n 

Consider, following [6], a two dimensional deformable medium. An explicit example 
could be a thin sheet of rubber or plasticine or any putty-like modeling material. Imagine 
that, at t — 0, an initial 'velocity' field (time rate of change of deformation) is specified for 
the medium. This will correspond to an expansion, a shear and a rotation (ESR) of the 
medium. Alternatively, one may specify (in place of initial velocity field) initial conditions 
on the ESR themselves. It may be recalled that the study of the evolution of the ESR 
variables of a medium is important in many areas (such as General Relativity). To proceed 
further, one needs to derive the evolution equations of the deformation in terms of the 
ESR variables. This is discussed below. Let us denote the deformation of the medium in 
terms of the deformation vector £' (where i = 1,2), which represents the vector joining any 
two infinitesimally separated points of the medium. The time rate of change of £ l may be 
expressed, for small time intervals, as (see Appendix A) 




(2.1) 
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where B l j(t) is an arbitrary second rank tensor characterising the time evolution of the 
deformation vector. Differentiating (12.11) with respect to time and using (12 .ip to eliminate 
£*, we obtain 

(dB^ \ 
-g> + B\B k A e = e- (2.2) 

Before we write down the evolution equations for the ESR, let us first define mathematically 
the quantities of interest. The arbitrary second rank tensor B l j can be decomposed into its 
trace, symmetric traceless and antisymmetric parts which will correspond to the isotropic 
expansion (trace scalar), shear (symmetric traceless tensor) and rotation (antisymmetric 
tensor). This decomposition is given as: 

Bij = 2®^ij + a ij + U V> ( 2- ^) 

where 9, and LUy represent, respectively, the expansion, shear and rotation. We can 
explicitly write (12. 3p in terms of the following 2x2 matrices: 



1 / \9 , 

:95 u = 2 , a. 



( 



2 \9 



ij 



2° ) \ 




uo tJ =\ ^ , (2.4) 



a + a x 
cr x -cr 4 

where we have denoted the shear components as a + and cr x , and the only rotation component 
as u. Thus, these four quantities 9, a + , a x and u characterise any deformation of the two 
dimensional deformable medium as a function of time. One can find useful discussions on 



the geometrical interpretation of each of these quantities in 
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B. The evolution (Raychaudhuri) equations 

We now use the inputs from the previous subsection to rewrite the evolution equation for 
B^ as four coupled, nonlinear, first order equations involving the dependent variables 9, a + , 
<r x and uj. Before that, of course, we need to write down a general expression for £\ This, 
we assume to be (in the linear approximation) of the form 

e =-K i j e-vi\ (2.5) 

where Kij and (3 represent the stiffness and viscous damping in the medium, respectively, 
with of the form 

(k ~\~ k-i. k x \ 
(2.6) 
k x k — k + J 
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Using ([23]) in (Q leads to 

B) + B\B) + K) + PB) = 0. (2.7) 

Finally, using the definition f)2.3p in (12. 7p . one obtains the differential evolution (in time) of 
the ESR as given by the equations: 

Q + ~0 2 + /39 + 2(a 2 + + a 2 x -uj 2 ) + 2k = 0, (2.8) 

a+ + (J3 + 6) a + + k + = 0, (2.9) 

& x + (J3 + e)a x + k x = 0, (2.10) 

w + (/3 + ^)cj = 0. (2.11) 

The equation for 6 is known as a Ricatti differential equation in the mathematics literature. 
It is nonlinear, first order and its solutions are known. Note that the equations for cr + , a x 
and uj are structurally similar. We now look for analytical solutions of these equations and 
try to identify generic features which appear in the solutions. 

C. Analytical solutions 

To find solutions, it is useful to look at specific cases before we attempt the most general 
situation. To this end, we consider the following cases : 

1. Kij = and (3 = 0, 

2. = k5ij and (3 = 0, 

3. = k5ij and (3 ^ 0. 

The equations f l2.8l) -( [2~TTi) for each of the above cases can be obtained easily, and analytical 
solutions can then be derived. We now discuss below each case in detail. 

Case 1: K {j = and (3 = 



We first consider the following pair of equations, 

9 + ^9 2 + 2I = 0, 
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(2.12) 



I + 291 = 0, 



(2.13) 



where / = o\ + cr x — u 2 . It is clear from the equation for / that its general solution would be 
of the form / = I e~ 2 f e ^ dt . This means that if I — at t — 0, it would remain zero forever. 
Further, I can not change sign during its evolution. The scheme for solving equations (I2.12p 
and ( 12.131) would be to replace 9 in (12.121) by — |j, then solve it for / and use it to find 
9 from (12.131) . Once 9 is known, the equations for a + , a x and uj can be solved easily. It 
is therefore obvious that the solutions to the equations (12.121) and (12.1 3D depend on the 
conditions imposed on I. We use this fact to obtain the different types of solutions. The 
exact solutions for 9 , <r + , a x and u with different conditions on / are given in Table [H 



I 


e 


{ <T+, cr x , U } 


I > 


c(d+^) 




{ E, F, G} 


2 (D+«)2_i6 




(D+^)2-16 




/ = 




{ ff+Oi O"x0, ^0 } 


(i+V) 




/ < 


C(D+S1) 




{ E, F, G } 


2 (D+^) 2 +16' 









TABLE I: The analytical solutions for Case 1 in two dimensions. 



In Tabled! C, D, E, F and G are the integration constants, andwill be used throughout 
hereafter in the text to represent the integration constants for the analytical solutions of 
other cases as well. Further, we have used { } brackets in Tabled to express the solutions in 
a compact way for certain ESR variables which have a common factor. This notation has 
been used throughout in the text. For the present case, the integration constants can be 
defined in terms of the initial conditions on the ESR variables at t — using D = where 
Io = er+o + o"xo — (the extra subscript indicates the initial value of the corresponding 
ESR variable). The other integration constants can then be easily written as follows : 

1. For / > 0: 



C 



29 



D 2 -16 



{E,F,G} = {a +0 , <j x0 ,u }{D 2 -16). 



(2-14) 
(2.15) 
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FIG. 1: Regions of initial conditions for singular and non-singular solutions in two dimensions. 
2. For I < 0: 



C 



29n 



'D 2 + 16 s 



{E,F,G} = { a +0 , a x0 , co } (D 2 + 16), 



(2.16) 
(2.17) 



where D is redefined as J_ Io . For 1 = 0, the solutions in Table Hare given in terms of the 
initial conditions on the kinematical quantities themselves. 

We note from Tabled that for / > 0, the solutions will tend to a singularity in finite time 
whenever D < 4, i.e., 9q — 4/ < 0, while for I = 0, we have a singularity for 9 < 0. On 
the other hand, there is no singularity for l > 0. This dependence of the occurrence of a 
singularity on the initial conditions is summarized on the Iq-9q plane in Fig. [D In Figs. [2] 
and El the evolution of the ESR variables is shown for different initial conditions. The initial 
conditions given in Fig.[2]^a) satisfy Iq > 0, while for the initial condition of Fig. [2(b), Iq < 0. 
It may be noted that the two initial conditions differ only in uj q . In the second case (i.e., 
Iq < 0), the singularity condition 9 — 4/ < can never be satisfied. The negative sign of 
Iq can only occur if u is large enough (i.e., Uq > cr^ + cr^o)- Thus, it may be inferred that 
there is a critical value of ujq that makes the solution non-singular. In other words, one can 
say that presence of a minimum rotation can prevent the appearance of a singularity. This is 
explained further in Section Hi Dl In case the amount of rotation present is not large enough 
(i.e. Iq > 0), the solution can still avoid a singularity if the expansion is large enough to 
prevent a collapse. This minimum value of expansion variable is obtained again from the 




FIG. 2: Case 1 in two dimensions: (a) for Iq > 0, and (b) for Iq < 0. 

singularity condition. Figure [3] shows the evolution of the ESR variables for two values of 
8 (when 7 > 0). 
Case 2: Kij = kSij and (3 = 
In this case the equation fl 2 . 1 3 [) remains same, while f)2.12p has an additional term 2k on the 
left hand side. Therefore, the form of solution for I will be same as in the Case 1. The exact 
solutions for ESR variables with different possible conditions are then given in the Table [Til 



/ 


8 


{ CT+, CJ X , U) } 


I > 

a 2 = 16 + b 2 , p 2 = 16k 


pa cos(D+t£) 


{ E, F, G} 


2[a sin(£>+f )+&] 


[a sin(D+f )+6] 


1 = 


2 v / A : tan[v / fc(C - t)\ 


{ D, E, F } sec 2 [Vk{C -t)} 


I < 

a 2 = 16 + b 2 , p 2 = 16k 


Q pa sin(D+^) 


{ E, F, G } 


~ 2[acos(D+f )+&] 


[a cos(D+f )+b] 



TABLE II: The analytical solutions for Case 2 in two dimensions. 
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FIG. 3: Case 1 in two dimensions for two different initial values of expansion scalar t9 . 

The constants of integration in terms of initial conditions are given as follows 
I. For I > 0: 

b = 



D 

{E, F, G} 



tan- ( P - 



291 

Ai ' 2/ () , 



29 J ' 
{ cr+o, &xo, } (a sin-D + 6). 



(2.18) 

(2.19) 
(2.20) 



2. For 7 = 0: 



3. For I < 0: 



C 

{A E, F} 



tan 



-l 



(9, 



V2v / ^/' 
{ °"+0) cr xo, } 



sec 



D = tan 



29 



{E, F, G} = {(T +0 , a x0 , uj } (a cos D + b), 



(2.21) 
(2.22) 



(2.23) 
(2.24) 
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where b is same as given in equation (12.181) for / > 0. In this case, it can be easily noticed 
that all solutions become singular in finite time. 



Case 3: = k5ij and (3 ^ 

For this general case, the equations (12.1 2p and (12.131) are modified as follows 

9 + -9 2 + (39 + 21 + 2k = 0, 



^2.25) 



I + 2((3 + 9)1 = 0. (2.26) 

The equation (I2.26P has its general solution of the form Further, we 

have various special cases resulting due to different possible conditions on the combination 
of stiffness and damping. The exact solutions for the ESR variables are calculated and given 
in the Table IIHI The integration constants used in Table II III may then be defined in the 
following form 



1. For I > 0: 



(a) (3 2 > 4k 



b 2 < 16, 



, -fo 2I ( 



D = tanh" 1 



2(J3 + 9 Q ) 

{E,F,G} = {o+o, cr x0 , u; } (a sinhD - 6). 



(2.27) 



(2.28) 
^2.29) 



b 2 = 16, 



D = In 



2(P + 9 



{F,F,G} = { a +0 , a x0 , lo }[1- exp(C)] 1 



(2.30) 
(2.31) 



b 2 > 16 



D = tanh 



-i 



2(P + 9 



_p + by/l _ 

{E, F, G} = {cr+o, cr x o, ^o} (a cosh£> - b). 



(2.32) 
(2.33) 



10 



/ > 

(i) P 2 > 4k 
• b 2 < 16 

„2 _ 



16 - 6 2 , p 2 = 4(p 2 - 4fc) 
6 2 = 16 



b 2 > 16, a 2 = b 2 - 16 



(it) /3 2 = 4A; 



(m) /3 2 < 4fc 



a 2 = 16 + 6 2 , p 2 = 4(4A;-/3 2 



pacosh(D+^) „ 
2[asinh(D+f )-f>] ~~ ^ 



pexp(J>+f ) 
2[cxp(D+f )-6] 



pa sinh(D+^) 
2 [a cosh(D+f )-&] 



2[(f + ^) 2 -16] ^ 



pa cos(D+^) ^ 



2 [a sin(D+f )+&] 



{ <7+, cr x , } 



{ E, F, G} 
[asinh(D+f )-6] 



Lg, F, G} 
[ cxp (D+f)-b] 



{ §, F, G} 
[a cosh(D+f )-b] 



{ E F, G} 
I(f + ^) 2 -16] 



{ E, F, G} 
[a sin(D+f j+b] 



1 = 

(i) P 2 > 4k, a 2 = I3 2 - 4k 

(ii) P 2 = 4k 

(iii) p 2 < 4k, a 2 =4k- P 2 



a [cxp(af )+cxp 2C ] a 
cxp(at)-cxp(2C) ~~ <° 



2 P 



atanf(C-i) 



{ D, E, F}cxp(at) 
[cxp(at)-cxp(2C*)] 2 



{ D, E, F} 

(t-cr 



{D, E, F } sec 2 §(C-i) 



J < 

(i) P 2 > 4k 

a 2 = 16 + b 2 ,p 2 = 4(P 2 -4k) 

(ii) P 2 = 4k 

(Hi) p 2 < 4k 



pasinh(D+^) „ 
2[acosh(D+f )+b] ~ P 



C(f + ^) 
2[(f + ^) 2 +16] 



pa cos(D+f ) 
2 [a sin(D+f )+&] 



{ E, F, G} 
[acosh(D+f j+b] 



{ E F, G} 
I(f + ^) 2 +16] 



[a sin(D+f )+&] 



.2 _ l2 



6 2 - 16 , p 2 = 4(4fc - p 2 



TABLE III: The analytical solutions for Case 3 in two dimensions. 
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(b) (3 2 = Ak 



(c) (3 2 < Ak 



2. For 1 = 0: 

(a) (3 2 > Ak 



(b) p 2 = Ak 



(c) (3 2 < Ak 



3. For I < 0: 

(a) /3 2 > 4A; 



4(/3 + O ) 



C = (/3 + o ) 
{£, F, G} = {cr+o, o- x0 , uj } 



(D 2 - 64) 



(D 2 - 64) 



6 = 



D 

{E, F, G} 



tan 



Jo ( 2(P + 9 ) 2 p 2 " 
r \ h 27 0/ 



-i 



2(/3 + o ). 
{o- +0 , cr x0 , uj } (a sinL> + b). 



exp(2C) 



/3 + #o - a 



f3 + 9 + a 

{D,E,F} = { a +0 , a x0 , w } [ 1 - exp(2C) f. 



C 



{D,E,F} = {a +0 , a x0 ,uj }C 2 . 



C 

{D, E,F} 



2 - tan- ?!L] , 
a \ a J 

{ °"+0; CxO, } 

sec 2 <f 



b = 



Iof_ 8 + 2([3 + 9 ) 2 p 2 



~ 2(/5 + ^ ) 



D = tanlT 1 

p- by f -I Q _ 

{E,F,G} = {cr+o, a x0 , } (a coshD + 6). 
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(b) (3 2 = 4k, 



D 



4(/3 + #o) 



(2.49) 



C = (P + d ) 



(D 2 + 64) 
D 



(2.50) 



{ E, F, G } = { a+o, ct x0 , uj } 



(D 2 + 64) 
4 



(2.51) 



(c) /3 2 < 4k 




(2.52) 



(2.53) 



{E, F, G} = {a+o, cr x0 , uj }(a sin D + b). 



(2.54) 



In this case, the regions of the ESR variable space of initial conditions that lead to solution 
singularity in finite time is more complex. We identify first some classifying features of the 
solutions presented in Table [III For Iq > and (3 2 < 4k, the solution is always singular, 
while for (3 2 > 4k the solution is non-oscillatory. These non-oscillatory solutions may or may 
not be singular depending on further conditions. For example, for (5 2 > 4k and b 2 < 16, 
we have singularity when sinhD < b/a. For I = and f3 2 > 4k, the solution has finite 
time singularity when (3 + #o + a < 0, as can be easily checked. The remaining sub-cases in 
Table ITO1 can be analysed similarly. It may be noted that, for 1$ < 0, allsolutions are non- 
singular, and we have an oscillatory solution only when f3 2 < 4k. The evolution of the ESR 
variables have been computed from the analytical solutions and presented in Figs. H] and O 
For Iq > 0, the effect of damping in avoiding solution singularity is observed in Fig. HI The 
existence of singularity is evident from Fig. H^a) for f3 2 < 4k, while there is no singularity 
in Fig. H](b) for the initial conditions which satisfy (3 2 > 4k. It is interesting to note that 
for f3 2 > 4k, the expansion variable settles to a negative value implying that the medium 
contracts indefinitely, however without the appearance of a singularity. On the other hand, 
for Iq < 0, the presence of stiffness produces oscillations as seen in Fig. [5j However, there is 
a minimum value of stiffness which induces oscillatory behaviour and is given by k > (3 2 /4. 
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FIG. 4: Case 3 in two dimensions: role of damping in suppression of singularity. 




FIG. 5: Case 3 in two dimensions: effect of stiffness in generation of oscillations. 
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FIG. 6: Sequence of deformed configurations of a square element in two dimensions for initial 
conditions given in Fig. [2(a). 

D. Visualisation of deformations of the medium 

The deformations of a medium in two dimensions may be easily visualised by considering 
the deformations of, say, a square element defined by four points of the medium. For such 
an element, we compute the deformation vectors from the central geodesic (taken as the 
centroid of the square) to the four corner points of the square at successive time steps. 
Thus, the four points define the shape of the element at any time instant. 

Let us consider the first case, with initial conditions as given in Fig. 12(a)- The evolution 
of ESR variables in Fig [2(a) clearly shows a singular behaviour in finite time. In order to 
understand this behaviour in terms of the deformations of the medium, we have plotted the 
successive shapes of an initially square element in Fig. [61 It is observed that, as the element 
shrinks, its area approaches zero. However, when I > 0, the rate of change of area A ^ 
as A — > 0. The expansion variable in two dimension has the interpretation of area strain 
rate, i.e., 9 = A/A (see {(J). Hence 9 = A/A — ► — oo, signifying a singularity. This can be 
easily seen as follows. Using this definition in (12 . 8 [) - ( T2TTTT) (with Kij = and (5 = 0), one 
can express A = J~4Jq/~A^ + PA, where A is the initial area of the element, and P is a 
constant. In the limit A — > 0, A — > (2/A )y/T^. For the initial conditions of Fig. [2(b), the 
successive shapes of the square element is shown in Fig. It may be noted in this figure 
that the square initially shrinks when 9 < 0, but later expands when 9 > 0. The presence of 
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FIG. 7: Sequence of deformed configurations of a square element in two dimensions for initial 
conditions given in Fig. [2jb). 





(a) 



(b) 



FIG. 8: Schematic representation of deformation of a square element in two dimensions for negative 
expansion (a) without rotation, and (b) with rotation. 

sufficient rotation helps in avoiding the singularity situation, as explained through Fig. [8] (for 
no shear). In the case when the expansion is negative, and shear and expansion variables 
are zero, the element collapses isotropically, as shown in Fig. [S](a). However, when a slight 
rotation is added to the collapsing element, singularity is avoided, as observed in Fig. Mb)- 
It may be noted in Fig [8](b) that when a corner point of the square (marked by a circle) is 
closest to the central geodesic, rapid rotation of the element is observed. In case shear is also 
present, there is a minimum initial rotation required to avoid a singularity. This minimum 
value of initial rotation is that which makes Iq < 0. Finally, one may observe that, due 
to the positive tail of the expansion variable 9 in Fig. EJ^b), the element in Fig. [7] expands 
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FIG. 9: Sequence of deformed configurations of a square element in two dimensions for initial 
conditions given in Fig. [5jb). 



even as shear and rotation of the element cease. For the case of initial conditions shown 
in Fig. Mb), where the response of all the variables is oscillatory, the deformations of the 
medium is visualised in Fig. [91 It is interesting to note that, due to the presence of damping 
in this case, the element progressively shrinks in size. Further, periodic occurrence of rapid 
shear and rotation about the points, where the expansion variable changes sign is observed. 
In case the damping is absent, no shrinking of the area will take place. 



III. DEFORMATIONS IN THREE DIMENSIONAL MEDIA 



A. Expansion, rotation and shear 

The basic scheme of analysing the deformations in two dimensions extend to three dimen- 
sions easily. However, we now have nine independent variables here (one for expansion, five 
for shear and three for rotation). So, the level of complicacy is certainly higher as compared 
to the two dimensional analysis in the previous section. 

As before, we begin by writing down the general expression for B^: 

Bij — Q^^ij + a ij + ^ij, (3-1) 
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where the expansion (9), shear (cr^-) and rotation {uJij) are given, 



/ i 



-95. 



9 



o \e o 



\o o \9) 



011 012 013 
""12 0"22 023 
\o"l3 °"23 — (fH + <7 2 2) 7 



-UJ 3 UJ 2 



UJi 



uj 3 — uji 
y -cj 2 o;i 



(3.2) 



B. The evolution(Raychaudhuri) equations 

The evolution equations (i.e. Raychaudhuri equations) for the expansion, shear and rotation 
can be then obtained in three dimensions in flat space. The nine coupled nonlinear equations 
for 9 (expansion), <7 n , cr 12 , cr 22 , o 13 , o 2 3 (shear) and u 12 , u>i 3 , u 23 (rotation) are derived (for 
the case when Kij = k5ij and (3 ^ 0) in the following form, 

9 + -9 2 + 09 + 2(<7 n + a\ 2 + a 2 2 + a 2 3 + a 2 3 + a n a 22 - uf - u% - uf) + 3k = 0, (3.3) 
2 1 2 

0"ll + (/? + g0) °"H + g(^ll + °"12 + °"13 - w| - W3) - -(of 2 + ^23 + ^11^22 ~ W?) = 0, (3.4) 

2 1 2 

022 + (/? + o 6 *) °"22 + ,(<t|i + 0-22 + ^23 ~ ~ ^l) ~ oKl + °"l3 + °"ll°"22 - ^2) = °, ( 3 - 5 ) 



012 + (/? + -6* + an + 022) 0-12 + 013023 + ^1^2 = 0, (3.6) 
o 

2 

013 + (/3 + -9 - cr 22 ) 0-13 + 012023 + CJ1W3 = 0, (3.7) 
2 

023 + (/? + -6 1 - (Tii) (T23 + CT13012 + CJ2W3 = 0, (3.8) 

2 

u 1 + (f3 + -9 - on) wi - cj 2 012 - ^3013 = 0, (3.9) 
2 

^2 + (/9 + -9 - o 22 ) cj 2 - cjiOi 2 - u) 3 a 23 = 0, (3.10) 
2 

^3 + (/? + -6* + On + o 22 ) ^ 3 - ^1013 - ^2023 = 0. (3.11) 



It is possible to note some structural similar between subsets of the above set of equations. 
The equations for on and o 2 2 are similar and we can obtain one from the other by inter- 
changing and replacing indices appropriately. In the same way, the three equations for the 
off-diagonal components of the shear, i.e. o i2 , Oi 3 and o 2 3 are also similar in structure. A 
structural similarity can also be noted for the three equations representing the evolution of 
the components of the rotation. 
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C. Analytical solutions: some special cases 

In the following, we solve the above equations under the following assumptions: 013 = 
(723 = 0, ui = uj 2 = and u 3 = iv. We now use the definition J = a\ 2 — uj 2 . It may be 
checked that these assumptions are consistent, i.e., they remain satisfied for all times if they 
are satisfied at t — 0. Now we have only five variables representing the ESR, namely, 9, an, 
012, 022 and iv. The evolution equations for these along with the procedure followed to solve 
the equations are given in the Appendix B. We have obtained analytical solutions for the 
following two cases: 

1. For Kij = and (3 = 0, 

2. For K i:i = k5ij and (3 = 0. 
These are now discussed in detail. 

Case 1: K {j = and f3 = 



1. J > 0: 

The exact solutions for 6, an, a 12 , a 22 and uj are given below, 



e 



+ 



1 



2E (D + Et) 

[(D + Ety-i] 1 (t-cy 



{a n , 022 } 
W12, ^} 



e 



{F,-F} 



3 (t-Cj [(D + Et) 2 — 1] 
{G, H} 



[(D + Ety-iy 

where C and D are given in terms of initial conditions as follows, 

3 



C 
D 



3 [(011)0 + (022)0] - 0o ' 
2 6> + 3(0n)o + 3(022)0 



(3.12) 
(3.13) 
(3.14) 



3^(011)0- (022)0 ] 2 + 4 Jo' 
The integration constants E, F, G and H can be written in terms of C and D as 



(3.15) 
(3.16) 



E 



(i + e c)(D 2 -i) 



2CD 

F = [(011)0- (022)0 ](D 2 -l), 
{G,H} = {(a 12 ) ,cu }(D 2 -l). 



(3.17) 
(3.18) 
(3.19) 
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These six integration constants (i.e. C-H) satisfy the condition 

2(D 2 - 1) 2 J + F 2 — AE 2 = 0. 



2. J = 0: 



1 1 1 

+ T. ^7 + 



(t-C) (t-D) (t-E) 



0u 



022 



1 



1 



(t -D) (t- C) (t - E) 
2 1 1 



{Cl2, ^} 



(t -E) (t- C) (t - D) 
{F,G} 



(3.20) 

(3.21) 
(3.22) 
(3.23) 
(3.24) 



(t-D) (t-E)' 

where D, E, F and G are given below in terms of initial conditions while C is same 
as in (13.151) . 

3 



D 

E 



3(oii) + <V 
3 



3(022)0 + 0o 
{F,G} = {(a 12 ) , u }DE. 



3. For J < 0: 



(a) [( ( 7ii)o-(a 22 ) ] 2 <4|J | 



(3.25) 

(3.26) 
(3.27) 





{ an, 022 } 

{012, ^} 



2E (D + Et) 1 
l(D + Et) 2 + l] + Jt^C) 







I + {F,-F} 



3 (t-C) ' \(D + Et) 2 + l] ' 



(3.28) 
(3.29) 
(3.30) 



[{D + Etf + lY 

where C is same as in ( 13. 15ft and D are given in terms of initial conditions as 
follows 

20 o + 3(0-ii)o + 3(0-22)0 



D 



E 
F 



3 V[ (011)0 -(022)0 ] 2 " 4 Jo 
(1 + 6 C)(D 2 + 1) 



2CD 

[(012)0 -(022)0] (D 2 + l), 
{G,H} = {(0 12 ) O) w }P 2 + 1). 



(3.31) 

(3.32) 
(3.33) 
(3.34) 
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These integration constants satisfy the following condition, 

2(D 2 + l) 2 J + F 2 - AE 2 = 0. (3.35) 

(b) [( ( Tii)o-(^ 2 )o] 2 = 4|J | 



' + TT^ (3-36) 



(t-C) (t-DY 



fan, 022} = - 



0_ 1_ {E, —E } 

3 (t - C) (t - D) 2 



(3.37) 



{012,^} = ^^, (3.38) 

where C is same as in (13.151) while other integration constants are defined as 
below, 

2C 

E = {(ffn)o- (022)0 }D 2 , (3.40) 

{F, G} = {(<7 la ) , ^o}^ 2 - (3.41) 



(c) [((7u)o- (o 22 )o ] 2 >4|J | 

The solutions remain same as for the case of J > with the integration constant 
D given by (13.3ip . 

It may be noted that the regions in the space of initial conditions that lead to finite time 
singularity of the solutions are more complex in the three dimensional case. For example, for 
Jo > 0, the solution has singularity when either C > (defined in (13.151) ) or D < 1 (defined 
in (I3.16P ). On the other hand, for J < 0, we have singularity only when C > (defined 
in (I3.15P ). It may be observed from (13 . 1 5[) that C = —1/B 33 . One can interpret B 33 as 
the rate of expansion/collapse (depending on whether .B33 > or ^33 < 0) of an element 
perpendicular to the plane of deformation. Thus, when C > 0, the element collapses in 
the direction perpendicular to the plane of deformation. For the case with initial conditions 
shown in Fig. [TUT a). Jo > 0, C < and D > 1. These conditions imply a non-singular 
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9 = 0.2, («„)„- 0.08, (o 12 ) = 0.03, (o^ = - 0.05, % = 




e Q = -0.01, ( Ol1 ) = -0.06, (a 12 ) = 0.01, (a 22 ) Q = - 0.05, a> Q = 0.02 



0.3 
0.2 
0.1 



-0.1 
-0.2 



I I I I I I 


(b) 






\ " ' / / 








1 1 1 1 1 1 





10 20 30 



40 50 60 70 80 90 100 

t 



FIG. 10: Case 1 in three dimensions: similarity of solutions as also observed in two dimensions. 
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0.5 r 
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t 

8 = 0.25, («„)„- 0.04, (o la ) — 0.2, (c 22 ) = 0.03, ^ = 0.25 




FIG. 11: Case 1 in three dimensions: effect of rotation in avoiding singularity. 
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solution as discussed above. On the other hand, for initial conditions satisfying J < and 
C < 0, the results are shown in Fig. flDTb). As expected, the solution in this case also is 
singularity free. 

The effect of rotation in avoiding singularity has been discussed in detail for the two 
dimensional case before. In the case of three dimensional deformations also, rotation may 
play a role in rendering the solution non-singular. Note that C in (13.15P is unaffected by uj. 
Thus, the necessary condition for rotation to make the solution non-singular is that C < 0. 
Now, for J > 0, the solution is singular if D < 1. Thus, if Uq is chosen such that either it 
makes D > 1, or makes Jo < 0, singularity is removed. For a set of two initial conditions 
which differ only in the value of lvq, the results are shown in Fig. [TTJ It may be noted that 
here the singularity is removed by increasing ujq in Fig. ITTfb) to make Jq < 0. 



Case 2: K 



ij — k5ij and (3 = 



1. For J > 0: 

The exact solutions for 9, cr n , <t 12 , <j 2 2 and to are given by 



[(q + pcos2(D - Vkt)] 



{(Tn, 0"22} 
{0"12, ^} 



- - Vktaii(C -Vkt) + 
3 

{F,G} 



{E, -E} 



[{q + pcos2{D - Vkt)}_ 



(3.42) 
(3.43) 
(3.44) 



[(q + pcos2(D - Vkt)} 
where p and q satisfy p 2 — q 2 = 4k. The integration constants are given in terms of 
initial conditions as 



C = tan 



6, 



3Vk 



(cn)o + (^22)0 



In order to obtain other integration constants, we first determine q as 



1 



2Bn 



d2 

-4k + -±[6o-Vk tan(v / A?C) ] 2 + B { 



where B is defined as, 



Bn 



4k 



[(^ii)o- (^2 2 )o] 2 + 4J 



(3.45) 



(3.46) 



(3.47) 



The other constants can then derived as, 



D 



tan 



B (VktanC - 6 
2Vk(B - q) 



(3.48) 
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Now, E, F and G may be written in terms of C and D as follows 



E = (q + p cos 2D) [ (a n )o ~ (022)0 ] , 

{ F, G} = { (0-12)0, ^0} (q + P cos2£>). 
These integration constants satisfy the following condition, 

E 2 + 4 J [q + p cos 2D] 2 - 16A; 2 = 0. 



(3.49) 
(3.50) 

(3.51) 



2. For J = 



= Vk[tan(C - Vkt) + tan(£> - Vfc*) + tan(£ - >/fct)] , 

6 



&11 



O22 



{0-12, ^} 



+ Vktan(D -Vict) 

o 


— + Vktan(E - Vkt) 
3 

{F,G} 



(3.52) 
(3.53) 

(3.54) 

(3.55) 



cos(D -Vkt) cos(£ -Vkt) 
The integration constant C is same as for the case J > while others are given as, 



D = tan 



-1 



flp + 3 (aii)o 
3Vk 



E = tan -1 



6*0 + 3 (022) 
3Vk 



{ F, G } = { (0-12)0, u } cos D cos E. 



(3.56) 
(3.57) 



The analytical solutions for the present case for J < are also possible with different 
conditions. The solutions for the case when [ (on) — (o 22 )o ] 2 > 4| J | are same as for J > 0. 
It may be easily checked that all solutions with 7^ have finite time singularity. 

In the general three dimensional case, due to a large number of variables, it may not be 
always possible to obtain analytical solutions of the kinematical quantities. However, one 
can obtain the numerical solutions for such cases. 



IV. SUMMARY OF RESULTS 

In this article, our primary purpose has been to find exact, analytical solutions of the 
expansion, rotation and shear by solving the system of equations involving these variables. 
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We have also shown, in the cases of two and three dimensional media, how our solutions 
differ when the initial conditions on the variables change. 

A crucial result in both two and three dimensions, which, we believe, is a generic charac- 
teristic of the system of equations is the existence of a critical value for the initial condition 
on the rotation. Above and below this critical value we have observed the presence and 
absence of a finite time singularity in the solutions. However, it is also observed that in 
three dimensions a change in the initial rotation alone does not necessarily avoid a singu- 
larity. Additional necessary conditions are also required to hold good in order to avoid a 
singularity. 

Characteristic effects related to the presence of damping and elastic forces are also evident 
in the behaviour of the ESR. For instance, in the presence of damping alone the area of a 
typical element is found to shrink in time whereas for a purely elastic medium (without 
damping) the area is conserved over a period for oscillatory solution. To help the reader 
appreciate our results we have given a visual representation of our results by looking at the 
deformations of a two dimensional square element in time with typical initial conditions 
above and below the critical initial rotation. 

A notable fact about our analysis in two dimensions is the identification of specific regions 
in the space of initial conditions which give rise to either singular or nonsingular solutions of 
the ESR. We have not been able to do so in three dimensions because of a larger number of 
variables involved and also due to the absence of a unique variable T (as in two dimensions). 
It will definitely be more illuminating if we could devise a way of understanding the initial 
condition space in three dimensions in a manner similar to the two dimensional case. 

V. DIRECTIONS FOR FUTURE WORK 

Even though we have been able to solve the case of deformations in two dimensional media 
exactly, the more interesting and realistic case of three dimensional media in its full glory 
(i.e with all the nine variables) remains an unsolved problem. It is surely a daunting task 
analytically, considering the number of variables and the coupled and nonlinear nature of the 
equations. However, we hope to be able to extract some useful information by numerically 
solving these equations and making use of the analytical solutions described here as checks 
on our numerics. 
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A generalisation of the system of equations derived and solved in this article concerns 
the situation where the evolution of the deformation is not necessarily described by a linear 
equation. In other words, the equation for £ may, in general, contain nonlinear terms of the 
form (quadratic), £ l £ J '£ fc (cubic) etc. Such terms would however lead to the presence 
of £ in the equations for the ESR variables. A probable way to solve such a scenario is 
to include the equation for £ in the system of equations and solve for the whole set either 
numerically or analytically. It is likely that an analytical evaluation will not be possible 
unless we consider a largely simplified model. 

We have mentioned in the Introduction that such equations could describe a variety 
of mathematically similar situations. The next obvious step in our study of such evolution 
equations for flows, is, inevitably, the investigation of fluid flows in two and three dimensions. 
Attempts towards analysing fluid flows are on and will be described in future articles. 

It may be apparent to the reader that the sequence of investigations being carried out (i.e. 
for deformable media, fluid flow etc.) seem to provide an indication towards the analysis of 
the full set of Raychaudhuri equations in a curved spacetime geometry for geodesic and non- 
geodesic flows. That is a truly formidable task, which, as of now, is not even a well-defined 
problem. As a preamble to such an exercise, it may be worthwhile to look at a curved space 
generalisation of the deformation of two and three dimensional media and then fluid flows. 

Finally, of course, we would like to see that our understanding of the kinematics of 
deformable media as described in this article actually does correspond to deformations in the 
real world. Our attempts at visualisation through simulations might actually be realisable 
in simplistic experiments. Note that one of the central results we have obtained concerns 
the critical value of the initial condition on the rotation, on either side of which we have 
nonsingular and singular solutions. This is also a feature which was noted in a cosmological 
context, i.e., rotation can prevent singularity formation in cosmology, though, spinning 
models (like the Godel universe) had other problems and, moreover, they did not correspond 
to a realistic cosmology. It would be instructive, if, in some way, we could show, even in a 
toy experiment, the correlation between rotation and singularity formation. 
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Appendix A 

Consider, at a time instant t, two infinitesimally separated points of a medium connected 
by the deformation vector £*, as shown in Fig. O At a later time instant t + At, the 
points move an infinitesimal distance along the respective instantaneous velocity directions, 
as shown. One can then write 

C(t + At) = f (t) - v\x j ,t)At + v i (x j + t) At 
or C = />>'/' 
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v l (x J , t) At 




FIG. 12: Schematic representation of evolution of deformation vector 



Appendix B 

The equations fl3.3l) - fl3.11l) . for case 2 of three dimensional deformable medium having 
Kij = k5ij and f3 = with cri 3 = <r 2 3 = 0, uj\ = u 2 = 0, = u, reduce to the following set 
of equations, 

+ z\ + z\ + (9 - z x - z 2 f + 2 J + 3k = 0, (B-l) 
Z! + 4 + J + k = 0, (B-2) 
z 2 + z\ + J + k = 0, (B-3) 
j + 2( Zl + z 2 )J = 0. (B-4) 



where zi = an + |# and z 2 = cr 22 + The equation (1B-4|) has the solution of the form 



J = J e- 2 f (zi+Z2)dt . Further, subtracting (jED and flB^l) from (jgTjl . one can calculate 
21 + ^2, and subtracting ( 1B-3I) from ( 1B-2I) gives 2:1 — 2:2, from which the different ESR variables 



can be calculated. In particular, for case 1 for three dimensional deformable medium, the 
equations (IB-lj) - (1B-3j) contain no term corresponding to the stiffness, and the solution for 



J remains same. 
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